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t— 5 This paper considers N x N matrices of the form Ay = A + jB, where 

A is self-adjoint, 7 G C and .£> is a non-self-adjoint perturbation of A. We 
obtain some monodromy-type results relating the spectral behaviour of such 
p 1 matrices in the two asymptotic regimes I7I — > 00 and I7I —> under certain 

C/3 assumptions on B. We also explain some properties of the spectrum of Ay 

for intermediate sized 7 by considering the limit N — > 00, concentrating on 
properties that have no self-adjoint analogue. A substantial number of the 
results extend to operators on infinite-dimensional Hilbert spaces. 
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1 Introduction 

Let A be a (possibly unbounded) self-adjoint operator acting in the Hilbert space 
H and let A y = A + 7I? where 7 G C and B is a bounded operator on 1-L. Many 
papers have been written about the spectral properties of the self-adjoint operators 
Ay when B = B* and 7 G R, the main techniques used including variational 
inequalities and perturbation expansions; see [5], [6J, [8] and many further references 
there. In this paper we concentrate on more general B and assume that 7 is 
complex. Our main concern is to describe phenomena that have no self-adjoint 
analogues, an issue that has been curiously neglected. A recent paper of Rana 
and Wojtylak, [7], is closer to this one, but there is little technical overlap. The 
interplay between the asymptotic regimes I7I — > and I7I — > 00 is a main focus of 
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interest, but we also explore some spectral phenomena that arise for intermediate 
values of 7. 

As well as being of intrinsic interest, operators of this type are relevant to non-self- 
adjoint Schrodinger and wave equations, for which the evolution is contractive as 
a function of time. In such situations every eigenvalue of A y lies in an appropri- 
ate half-plane and the eigenvalue determines the energy and rate of decay of the 
associated eigenstate of the system. From Section [6] onwards we study rank one 
perturbations. As well as providing a range of phenomena that must be included 
in a more general theory, this is of direct relevance to the study of non-self- adjoint 
boundary conditions for Schrodinger operators in one dimension. The relevant per- 
turbations of the Schrodinger operators are singular, but, if one considers instead 
the resolvent operators, the perturbations are rank one and bounded. 

General considerations from perturbation theory imply that the set 1Z of (7, A) G 
C 2 such that A is an isolated eigenvalue of Ay with finite algebraic multiplicity 
is a Riemann surface that may have branch points where the multiplicity of the 
eigenvalue is greater than 1; see [I]. If B is relatively compact with respect to A 
then 

Spec(A 7 ) = Ess(A) U {A : (7, A) G 71} 

for every 7 G C. Our goal in this paper is to understand how the geometrical 
structure of 7Z depends upon some simple generic assumptions about A and B. 

In much of the paper we assume that 7i has finite dimension N. We assume that 
A is self-adjoint and that B is sectorial. The coupling constant 7 is restricted by 
the requirement that Im((A 7 /, /)) > for all / G %] this is equivalent to assuming 
that iAj is dissipative in a standard sense; see [U Section 8.3]. Further assumptions 
on B are made as necessary. In the particular case B = B* > 0, which motivated 
our initial interest, we assume that < arg(7) < it. Theorem 26 and Example [29] 
show how a substantial part of the spectrum of a large matrix may sometimes be 
approximated by using a carefully chosen matrix that is much smaller. Section [7] 
focuses on spectral properties of v4 7 that are best understood by considering the 
limit iV — > 00. 



2 Sectorial operators 

The truncation of an operator A on H to a closed subspace K is defined by A^ = 
PAP\/c, where P is the orthogonal projection of % onto /C. We will need the 
following lemma. 

Lemma 1 If A = A* , al < A < bl and A^ denotes the truncation of A to K, then 
a/" < A* < blK 
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Proof We use variational methods. The hypotheses imply that 

a < M{(AfJ):fen and ||/|| = 1} 
< m£{(Af,f):feK and ||/|| = 1} 
= inf{(^/,/):/G/C and ||/|| = 1}. 

Therefore afl < A\ The other half of the proof is similar. □ 

A bounded operator D on the Hilbert space H is said to be sectorial if there exist 
'sectorial constants' o~\, a 2 such that —ir/2 < —o~% < < a 2 < 7r/2 and 

{(Df, /) : / e H} C {0} U {z : z ^ and - <j x < arg(^) < a 2 }. (1) 

The theory of sectorial operators has a long history; see Sections VI. 1.5 and VI. 3.1 
of The following lemma is adapted from [3J Theorem VI. 3. 2], but we include a 
proof for completeness. 

Lemma 2 If D is a bounded sectorial operator on % and f 6 % then the following 
are equivalent. 

(i) (Df, f) = 0; 

(ii) ((D + D*)f, f) = 0; 
(Hi) Df = 0; 
(iv) D*f = 0. 

If K, — Ker(D) then K. and K. 1 - are invariant under D and D* . Moreover D\^ = 
D*\k — 0. Both D\fc± and D*\^ are one-one with ranges that are dense in K, L . 
The truncation D^ of D to K 1 - may be written in the form 

D^ = X 1/2 {I* + iE)X 1/2 (2) 

where X is the truncation of (D + D*)/2 to K 1 - , 1^ is the identity operator on K, 1 - 
and E is a self-adjoint operator on K, 1 - satisfying 

- tan((Xi)^ < E < tan(a 2 )/\ (3) 

where o~\, o~i are the sectorial constants of D. 

Proof 

(i) implies (ii). This uses (D*f,f) = Jb~JJ). 

(ii) implies (hi) and (iv). We write D = D +iDi where D = D$ > and D\ = D\. 
The sectorial condition is equivalent to —kiD < D\ < k 2 D Q where k r = tan(oy) for 
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r = 1, 2. If (i) holds then (£>„/, /) = 0, so ||I>J /2 /ll 2 = (A)/, /) = 0. This implies 

1 /2 

that £) / = 0, and hence that D f = 0. Since < D\ + /cD < 2kD , we a ls° 
have (pi + jfcD )/, /) = 0, hence (A + kD^) 1 ' 2 f = and then (A + fc£> )/ = 0. 
Therefore A/ = 0. We conclude that Df = and D*/ = 0. 

(iii) and (iv) separately imply (i). Both are elementary. 

The property (iii) implies that D\jc = 0. The property (iv) together with the 
general identity 

Ran(L>) = (Ker(A)) 1 

implies that Ran(D) is dense in K, x . The corresponding statement for D* has a 
similar proof. 

We have, finally, to prove (|2| and O). Without loss of generality we assume 
that K, = and omit the symbol . The operator X = D is then one-one with 
dense range T> in The inequalities —kiD < D\ < k 2 D are equivalent to 
—k\I < E < k 2 I where E = D^^DiDq 1 ^ 2 is initially defined as a quadratic form 
on T>. This yields (J3j. The bounds on the form E imply that it is associated with 

1/2 1/2 J L 

a bounded linear operator on T-L. We then have D\ = D ED and hence (2). □ 



Corollary 3 If D is sectorial and S is bounded then the following are equivalent. 



(i) SDS* = 0; 

(ii) SD = 0; 

(iii) SD* = 0. 



Proof Assuming (i), (SDS*g,g) = for all g e U. Therefore (Df, f) = for 
all / e Ran(»S'*). Lemma [2] now implies that Df = for all / 6 Ran(5**). Hence 
DS* = and (3) holds. The proof that (i) implies (ii) is similar and the proofs 
that (ii) and (iii) separately imply (i) are elementary. □ 

The remainder of this section is of independent interest, but it is not used elsewhere. 
Given constants <J\, a 2 such that — 7r/2 < — <j\ < < ct 2 < n/2, the set of all 
bounded operators D on the Hilbert space H such that Q holds is a proper closed 
convex cone, which we denote by C ai}(T2 . We say that a non-zero operator C lies 
in dC aija2 if C = A + B and A, B 6 C aija2 imply that there exist non- negative 
constants a, (3 such that A = aC and B = (3C. The set of all positive multiples 
of such an operator C is called an extreme ray of C CTl)C r 2 - 

Lemma 4 Let A, B,C G C CTl)(J2 and C = A + B. Then Ker(C) C Ker(A). In 
particular rank(C) = 1 implies A = or rank(A) = 1. 
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Proof The assumptions imply that 



C + C* = (A+ A*) + (B + B*). 

and then 

< A + A* < C + C*. 

Therefore ((C + C*)f,f) = implies ((A + A*)f,f) = 0. Lemma [2] now implies 
that Ker(C) C Ker(A). □ 

Theorem 5 ^ Let C ai ^ 2 be the cone defined above. Then a non-zero operator 
A G dC aii(T2 if and only if Af = a(f,e)e for all f G %, where e G H satisfies 
\\e\\ 7^ and a = e~ lcri or a = e l<72 . 

Proof Given A G C ai>a2 , let /Ci = Ker(A). If KL± has dimension greater than 1, 
then by applying the spectral theorem to the self-adjoint operator E in ([2]), one 
may write /C^ = /C2 © /C3 where /C2 and /C3 are non-zero orthogonal subspaces that 
are invariant with respect to E. One then has a block decomposition 



P + iE 



h + iE 2 
h + iE 3 



in an obvious notation. Moreover I2 + 1E2 and J3 + iE% both lie in C ai ^ 2 with 
respect to the relevant Hilbert spaces. It follows that A = A 2 + A 3 where A 2 and 
A 3 have the following block decompositions with respect to H = JCi © /C 2 © £3- 


A 2 = A\' 2 \ I 2 + iE 2 ) A; /2 , 






A 3 = A 1 ^ (0 ] AT 

I 3 + iE 3 



1 /2 

The factors A do not change the sector in which the numerical range lies, so 
A 2 , A 3 G C ai>a3 and A ^ dC aur72 . 

Conversely if K,± is one-dimensional then A has rank 1 and it is of the form Af = 
{f, ^i)e 2 for some non-zero vectors e±, e 2 and all / G H. Since 

Ker(A) = {/:(/, e x ) = 0}, 
Ker(A*) = {/:</,e 2 > = 0}, 

Lemma |2] implies that Ker(A) = Ker(A*), from which one may deduce that e 2 is a 
multiple of ej. An easy calculation using Lemma |4] shows that A is in an extreme 
ray if and only if the argument of a has one of the two stated values. □ 
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3 Cyclicity 



This section generalizes the notion of cyclic vector to perturbations of an operator 
that have rank greater than 1. 

Theorem 6 Let A be a (possibly unbounded) self-adjoint operator acting in the 
Hilbert space H and let B, X be two bounded operators on H. Then the following 
conditions are equivalent. 

(i) Xe iAt B = for all t G R; 
(it) X(zl - A)- 1 B = for all z £ Spec(A); 

(in) Xe< A+ ^B = for some (equivalently all) 7 G C and all t G R. 
If A is bounded the above conditions are also equivalent to 

(iv) XA n B = for all n>0. 

Proof We use a number of standard theorems and formulae from the theory of 
one-parameter semigroups; see [Tj, Sections 8.2, 11.4]; in finite dimensions many 
of these can be derived more directly. We first observe that A = A* implies that 
there is a one-parameter group with generator iA; following the usual convention 
we write this in the form e lAt , where t G R. The boundedness of B implies that 
there is a one parameter group, which we denote by e^ A+lB ^, whose generator is 
A + ^B. 

(i)=^(ii). This follows directly from the following formulae, the integrals being 
convergent in the strong operator topology. If \m(z) < then 



-1 



= 1 




If Im(z) > then 




If z G R\Spec(A) and e > then 



(zl - A) 



-1 



lim(zl + iel — 



(ii)=^(i). This uses the formulae 



{si =f %A) 



-n-l 



(-l) n d 



(si =f iA) 



-1 



n! ds 
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The formulae are valid for all positive s, t and n and the limits may be taken in 
the strong operator topology. Both formula may be proved by using the spectral 
theorem, but they are also valid at the semigroup level. 
(i)=r*(iii). Assuming t > 0, this uses the formula 



e i(A +1 B)t = Q iAt + f e iA(t- S ) ijBe iAs dg 
Js=0 

t 

iA(t-8)' p>JA(s-u)- tjJ.Au 



+ / / e tA ^- s) i 7 Be lA(s - u) i 7 Be lAu duds + . . . , 

Js=0 Ju=0 

the integrals and series being convergent in the strong operator topology for all 
7 G C. The proof for t < is similar. 

(iii)=£-(i). If (iii) holds for some 7 G C then (i) follows by using the formula 

Q iAt _ e i(A+-yB)t _ f e i(A+ 1 B)(t-s) i ^ Be i(A+~/B)s ^ 
Js=0 

+ f f e i(A+ ^ )( *- s) z7 J Be 4 ( A+ ^( s - u ^ 7 Se i ( A+ ^)"d M d S + ... 

Js=0 Ju=0 



(i)<^>(iv). These use 



d n 

(iAT = W J M 



c 



t=o n=0 n - 



E 



both limits being in the operator norm. □ 

In the context of Theorem [6j we say that the bounded operator B is cyclic for A 
if the conditions of the following corollary hold. 

Corollary 7 Let A be a possibly unbounded self-adjoint operator acting in the 
Hilbert space H and let B be a bounded operator on H. Then the following condi- 
tions are equivalent. 

(i) Whenever any of the equivalent conditions of Theorem [#| holds for some 
bounded operator X on H, it follows that X = 0. 

(ii) If one defines 

C 2 = lm{e iM Bv : tell and v EH} 
then £2 is dense in %. 

(iii) If one defines 

£ 3 = \m{(sl - A)' l Bv : s <£ Spec(A) and v E U} 
then £ 3 is dense in %. 
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(iv) Assuming that A is bounded, if one defines 

£ 4 = \in{A n Bv : n = 0, 1, 2, . . . and v G H} 
then £4 is dense in %. 

Proof (i)=^(ii). If (ii) is false then the Hahn-Banach theorem implies that there 
exists a non-zero <fi G H such that (<p,v) = for all v G £2- If one defines 
Xv = (v,(fi)(f) for all v G W then one sees that Xe lAt Bv = for all v G % but 
X 7^ 0, so Theorem (6Fi) is false. 

(ii)=^(i). If Theorem[6ti) is false for some non-zero X G C(T-L) then £ 2 ^= Ker(X) 7^ 
so (ii) is false. 

The proofs that (i)-^(iii) and (i)-^(iv) are very similar. □ 

In the following theorem and elsewhere we use the notations C + = {z G C : 
Im(z) > 0} and C_ = {z G C : Im(^) < 0}. If B is a sectorial operator with 
sectorial constants a%, 02 we define 

S B = {0} U {7 G C : 7 ^ and a x < arg(7) < vr - a 2 }. (4) 

The condition 7 G Sb implies that "f(Bf, f) G C + U {0} for all / G %. 

Remark 8 The conditions in Corollary [7] only depend on B via the closure of its 
range TZq = {Bf : f G %}. In particular if A and B are both bounded, then B is 
cyclic for A if and only if the linear span of [J r>0 A r TZ is dense in %. Given m G N, 

let lZ m be the orthogonal complement of IJ^o 1 ^"T^o in [J^lo A r 7Zo. Then lZ m are 
orthogonal subspaces and B is cyclic for A if and only if the closure of the sum of 
{^m}m=o ^ s dense on %. One may use these subspaces to represent A as a self- 
adjoint block tridiagonal matrix. If B is sectorial and B r ^ s is its associated block 
matrix, then i?o,o is the truncation of B to TZq and all other entries B T ^ S vanish. If 
% is finite-dimensional, one only has a finite number of non-zero subspaces. □ 

Remark 9 The conditions in Corollary [7] are close to those used in the block 
Krylov subspace method of numerical analysis. Case 4 corresponds to the standard 
version of the method while Case 3 corresponds to the rational version. □ 

Theorem 10 Suppose that B is sectorial and that 7 G Sb- If B is cyclic for A and 
A is an eigenvalue of A^ then A G C+. If M = rank(I?) < 00 then the geometric 
multiplicity of X is at most M. 

Proof Suppose that 7^ / G Dom(v4 7 ) and Af + ^Bf = A/. By calculating the 
imaginary part of 

(AfJ) + ( 1 BfJ) = X(fJ) 
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one deduces that either A G C + or A G R and Im^jBf, /)) = 0. Since jB is 
sectorial it follows that (7!?/,/) = 0. Lemma [2] now implies that Bf = B* f = 0. 
Therefore Af = Xf and e iAt f = e iXt f for all t G R. Therefore B*e iAt f = for all 
t G R and 

(f, e- iAt Bv) = 

for all t G R and all v G "H. Since -B is cyclic for A it follows by Corollary [7^ii) 
that / = 0. The contradiction implies that A G C+. 

If (A + jB)f = Xf then (A/ - A)f = jBf. Since A G C+, A £ Spec(A) and 
/ = (A/ — A)~ 1/ -/Bf G (A/ — A)~ 1 B'H, which is a linear subspace with dimension 
at most M. □ 



4 The main theorems when N < oo 

In this section we suppose that iV = dim('H) < 00 and put M = rank(S) where 
B is sectorial. Our goal is to describe how the spectrum of A 1 = A + ^B depends 
on 7, assuming that 7 6 5b as defined in Q, and in particular the relationship 
between the spectral asymptotics for small and for large 7. 

Under the above assumptions it is elementary that lm(((A + t'yB)f, /)) is a mono- 
tonically increasing linear function of t G (0, 00), as is Im(tr((A + t7-B))). Combin- 
ing these observations with known variational results for B = B* > and 7 > 0, 
leads to the conjecture that the imaginary part of each eigenvalue of A + tjB also 
increases monotonically as a function of t. The following example demonstrates 



that this is false. It also illustrates the results in Theorem 19 Example 1.5.7 of 



PQ , which is even simpler, provided one of the motivations for the present study. 

Example 11 Let A be the 5x5 diagonal matrix with eigenvalues X r = r for 
1 < r < 5, and let A 7 = A + jB where B is the rank 2 operator 

Bf=(f,e 1 )e 1 + (f,e 2 )e 2 

for all / G C 5 , where e x = (2,2,2,2,2) and e 2 = (3,3,-2,-2,-2). Figure Q plots 
the eigenvalues of A 1 for 7 = te 10 , where < t < 00 and 9 = 37r/8. The eigenvalues 
converge to the eigenvalues of A as t — > 0. Two of the eigenvalue curves diverge as 
t — > 00, while the other three converge back to the real axis. □ 

We shall need the following conditions. Apart from (HI), each is generic in the 
sense that it holds for a dense open subset of operators of the relevant type. 

(HI) dim('H) < 00, A = A* and B is sectorial. 
(H2) The operator B is cyclic for the operator A. 
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Figure 1: Spectral curves described in Example 11 



(H3) All of the eigenvalues of A have algebraic multiplicity 1. 

(H4) All of the non-zero eigenvalues of B have algebraic multiplicity 1. 

(H5) All of the eigenvalues of the truncation of A to the kernel K. of B have 
algebraic multiplicity 1. 

Theorem 12 Let 7 e S B - If (HI) holds and Z = i(A + 7E) then \\e zt \\ < 1 for 
all t > 0. Given (HI), the condition (H2) holds if and only if there are constants 
M > 1 and c > such that 

\\e zt \\< Me~ ct (5) 

for all t > 0. Given (HI) and (H2), one can put M = 1 in if and only if 
Ker(5) = {0}. 

Proof It follows directly from (HI) that Z is dissipative for every 7 6 Sj and 
hence that e zt is a contraction semigroup for t > 0. If (H2) also holds then every 
eigenvalue A of Z satisfies Re(A) < 0, and an application of the Jordan form 
theorem yields □ 

Example 13 Suppose that A and B satisfy (Hl-5) and that every eigenvalue A 
of A satisfies A > 0. Define the operators A and B on H © H by A(f © g) — 



10 



(Af) © (cAg) and B(f © g) = (Af) © (cBg), where c> 0. Then A and B satisfy 
(H1-H5) for almost all such c, but not for c = 1. The proof uses Lemma [TJ □ 

We assume (HI), (H2) and that 7 G Sb throughout the section, so that we can use 
Theorem 10 We make constant use of the polynomial 

p(7, A) = det(A + 7_B - XI). (6) 

The 7-dependence of the spectrum of A^ depends on an analysis of the algebraic 
surface 

n = {( 1 ,\)eS B xC + :p( 1 ,X) = 0}. (7) 



We will use the following classical facts. 

Proposition 14 If X is an N x N matrix and q(X) = det(X — XI) then q is a 
polynomial with degree N and the following are equivalent. 

(i) Every eigenvalue of X has algebraic multiplicity 1; 

(ii) Every root X of q is simple; 

(Hi) There are no simultaneous solutions of q(X) = q'{X) = 0; 

(iv) The discriminant of q is non-zero. (The discriminant of a polynomial q is a 
certain multiple of the square of its Vandermonde determinant, and may be 
written as a homogeneous polynomial with degree 2N — 2 in the coefficients 
ofq-J 



Since the zeros of p(0, X) all lie on the real axis, the following lemma can often be 
used to reduce the determination of the zeros of ^(7, A) in to a lower dimensional 



problem. See Lemma 20 The right-hand side of <\8h, usually without the % is called 



the relative determinant of Ay and A. 



det((/ + 7 (A-A/)- 1 fi) h ) 



Lemma 15 One has 

KM) 

where ^ denotes the truncation of the operator to the range of B* 
Proof This is a combination of two identities 

= det(/ + 7 (X-A/)-> B ), 
= det ((J + 7(A - Xiy 1 B) ii ) . 
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The first equality is obtained by calculating the determinants of both sides of the 
identity 

A + 1 B - XI = (A - XI)(I + i(A - \iy 1 B). 

The second equality is proved by writing I + / y(A — XI)~ 1 B as a 2 x 2 block matrix 
using the orthogonal decomposition 

H = Ker(B) ©Ran(5*). 

□ 

Lemma 16 Given (H2) and (H3), there exists a finite set F\ C Sb, such that A^ 
has N distinct eigenvalues, each with algebraic multiplicity 1, for every 7 G Sb\F\. 
If{j,X)enandj^F 1 then^(j,\)^0. 

Proof The eigenvalues of Ay are the roots of the polynomial <7 7 (A) = p(j, A), which 
is of degree N in A with leading coefficient (— 1)^. The eigenvalues of v4 7 lie in C + 



by Theorem 10 They all have algebraic multiplicity 1 if and only if the discriminant 



of g 7 is non-zero, by Proposition 14 The coefficients of g 7 are polynomials in 7, 
so the discriminant is also a polynomial r in 7. The hypothesis (H3) implies that 
r(0) 7^ 0, so r is not identically zero, and it has only a finite number of roots. The 
first part of the proof is completed by putting F\ = {7 G Sb '■ r(j) = 0}. The 



proof of the final part of the theorem uses Proposition 14 again. □ 



Lemma 17 Given (H2) and (H4), there exists a finite set F2 C Sb, such that if 
(7, A) Ell and-i £ F 2 then ^(7, A) ^ 0. 

Proof One may evaluate p(j, A) by combining an orthonormal basis of Ker(B) 
with a set of M eigenvectors associated with the non-zero eigenvalues P%, ■ ■ ■ , (3m 
of B. If one does so then one sees that qx(j) = p(y, A) is a polynomial with degree 
(at most) M in 7 whose leading coefficient is det^ — Xfl) Yl r=1 (3 r , where A^ is 
the truncation of A to Ker(S) and P is the identity operator on this subspace. 
Since is self-adjoint and A G C + , the determinant is non-zero and the degree of 
q x is M. 



One may see as in the proof of Lemma 16 that the roots of q\ are all distinct if and 



only if a certain polynomial r(A) is non-zero. If p(y, A) = and ^(7, A) = then 
r(A) = 0. The set G\ of roots of r is finite provided r does not vanish identically. 
Assuming this, 

F 2 = {7 G C + : (7, X) Ell for some A G d} 

is also finite and §^(7, A) 7^ for all (7, X) en such that 7 ^ F 2 . 

The polynomial r is not identically zero provided the M solutions 7 of det(A + 
jB — XI) = are distinct for all large enough A G C + . This is true if and only 
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if the solutions s of \~~ N det(A + sXB — XI) = are distinct for all large enough 
A G C + . These solutions converge as |A| — > oo to the solutions of det(sB — I) = 0, 
which are Pi 1 ,..., (3 m- They are distinct by (H4). □ 



The following lemma will be used in the proof of Theorem 19 



Lemma 18 Let L be a bounded operator on/H — Hi® H2 with block matrix 

L-( P Q 
\R S 

where the entries satisfy \\Q\\ < c, \\R\\ < c, US' -1 !! < l/(2c) and ||-P _1 || < e < 
l/(2c). Then L is invertible and 



Proof If one puts X 



L- 1 



P 

s 





o s- 1 



and Y 



< 2e. 



Q 
R 



then ||Y|| < c and 



\X 1 || ^ V(2c)- Therefore \\YX 1 \\ < 1/2 and the perturbation expansion 

oo 

{x + y)- 1 = x~ l ^(-yr 1 )" 



n=0 



implies that L = X + Y is invertible with \\L 1 || < 1/c. Moreover 



L- 1 





o s- 1 



< \\(X + Y)- 1 -X^W+s 



< || - X^YX- 1 + (X^YX-^iYX 
-(X- l Y X- l )(Y X 



-1\2 



< 2\\X- 1 YX- 



P^QS- 1 
S^RP- 1 



+ e 



+ e 



< 2e. 



□ 

We use the above results to connect the spectrum of for large and small 7. Let 
Q denote the set of all continuously differentiable curves g : [0, 00) — > C such that 
g(0) = 0, g(t) G Sb for every t > and #'(£) does not vanish for any t G [0, 00). Let 
F = F\ U F2 where F\ is defined as in Lemma [16] and F2 is defined as in Lemma [17) 
Let ^0 denote the set of all curves g G Q such that g(0) = 0, G Sb\F for 
every £ > and lim^oo \g(t)\ = 00. In the next theorem, one can impose stronger 
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conditions on g (e.g. C°° or real analyticity) and obtain similarly strengthened 
conclusions on the eigenvalue curves A r . 

Our main theorem below is an example of monodromy in the sense that we prove 
that certain one-parameter curves that avoid a finite number of singularities may 
have different end points even if they have the same starting point, provided they 
take different routes around the singularities; the difference is measured by an 
element of a permutation group. 

Theorem 19 Given (H2-5), let g G Qq. Then there exist N curves A r G Q such 
that Spec(A g ( t )) = {Ai(t), . . . , A7v(i)} for all t G [0, oo). One can choose the order- 
ing of these so that A r (0) = a r for all r G {1, . . . , N}, where a r are the eigenvalues 
of A written in increasing order. Assuming this is done, there exists a g-dependent 
permutation n on {1, ... , N} such that 

lim basM = ft 

t^oo g(t) 

for 1 < r < M , where (3 r are the non-zero eigenvalues of B written in any fixed 
order, and 

lim A 7r (A/ +r ) (t) = 5 r 

t— >oo 

for 1 < r < N — M, where 5 r are the non-zero eigenvalues of the truncation A^ of 
A to Ker(I?) written in increasing order (the eigenvalues are distinct by (H5). If 
g G Go is a real analytic curve then so are all the curves A r . 

Proof If t > then g{t) F, so the eigenvalues Ai,...,Ajv of H g ^ all have 
algebraic multiplicity 1. Perturbation theory implies that each eigenvalue of A 1 is 
an analytic function of 7 if 7 = g{t). Therefore each eigenvalue \{t) of A g u) is a C 1 
function of t, or real-analytic if g is real analytic. These perturbation arguments 
imply all the statements of the theorem that relate to the limit t — > 0. 

We next observe that p(g(t), X(t)) = for all t > 0. Differentiating this with 
respect to t yields 

^(g(t),\(t))g'(t) + ^(g(t),X(t))\>(t) = 0. 



By applying Lemmas 16 and 17 we deduce that A'(t) is non-zero for every t > 0. 

In order to prove the remainder of the theorem we need only find the asymptotic 
forms of the eigenvalues of iJ 7 as I7I — > 00 and apply the results to 7 = g{t) as 
t — > 00. The spectrum of A^ is a set rather than an ordered sequence and there is 
no reason for any ordering of the eigenvalues of A-y for large 7 to be related to the 
ordering for 7 = 0. 
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We start by describing the large eigenvalues of Ay. For every r G {1,...,M}, 
perturbation theory and (H4) together imply that B+ , y~ 1 A has a simple eigenvalue 
of the form 

/ir = / 3 r + 7 -i(A/ r ,/;> + o( 7 - 2 ) 

as |7| — > 00, where f r is an eigenvector of B associated with the eigenvalue {3 r , f* 
is an eigenvector of B* associated with the eigenvalue f3 r and we normalize so that 
(f r , f*) — 1- This implies that Ay has a simple eigenvalue of the form 



\ r = 1 p r + (Af r j;) + o{ 1 - 1 ) 



for all r G {!,..., M}. 



We next use Lemma 18 to describe the small eigenvalues of Ay. If one defines 
Hi = Ran(£>) and Hi = Ker(B) then C N = Hi © H2 is an orthogonal direct sum 
by Lemma [2} One may write 

D / C + jB* E 

where Bfi is the truncation of B to "Hi, C is the truncation of A to "Hi and A is 
the truncation of A to H.2- We now add hi to both sides where the constant k is 
independent of 7 and large enough to ensure that ||(A + < l/(2c), where 

c = || £7 1| + 1 = || -£7*|| + 1. Using the fact that B^ is invertible on Hi, we observe 
that 

e= \\(C + 1 B* + kI)- 1 \\ =0(| 7 |- 1 ) 



as I7I — > 00. Lemma 18 now implies that A + ^B + kl is invertible and 

= Oi 



(A + 1 B + kiy l - ( 





(A + A;/)" 1 



(9) 



for all large enough |7|. Every eigenvalue S r of satisfies 

2 < 2c < \5 r + k\ < \\A*\\+k. 

Since A^ is self-adjoint, a perturbation argument applied to ^ implies that there 
is an eigenvalue /i r of A + 7S such that 

K^ + A^-i-^ + A;)- 1 ! =0(| 7 r 1 ) 
as | 7 | — > 00. By combining the last two equations we obtain 



\fJL r ~ 5 r \ = 0(\j\ 



-1> 



as | 7 | — > 00. Moreover, the perturbation argument proves that fi r has the same 
multiplicity 1 as 5 r for all r G {1, . . . , iV — M}. 
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We have now described N distinct simple eigenvalues of A 1 for all sufficiently large 
|7|. Since A 1 is an N x N matrix there are no other eigenvalues. □ 

Simple continuity arguments show that two homotopic curves gx, g^ G Qq give rise 
to the same permutation tt. The fact that non- homotopic curves may give rise to 



different permutations is demonstrated in Examples 37 and 38 



5 Localization 



In this section we describe a procedure for approximating the spectrum of = 
A + jB in a given region of C+. We assume that A is a (possibly unbounded) 
self-adjoint operator on H, that /C is an auxiliary Hilbert space, that B = CD and 
that C : K — > H, D : H — > K. are bounded operators. 

We first note that the Birman-Schwinger method does not depend on self-adjointness 
of the perturbation. Numerically, the method is most useful when the dimension 
of /C is much smaller than that of H, but one need not assume that either is 
finite-dimensional. 

Lemma 20 If A ^ Spec(A) and 7 7^ then X is an eigenvalue of A^ if and only 
if— 1/7 is an eigenvalue of 

m(A) = D(A - \iy l C G £(£). 

If K is finite- dimensional then A ^ Spec(A) is an eigenvalue of A 7 if and only if 
the jointly analytic function 

p(7,A) = det (J + 7771(A)) 

vanishes. 



Proof We start with the identity 

A + 7CD - XI = (/ + ^CD{A - \iy l )(A - XI), 

both sides being regarded as linear maps from Dohi(t4) to H. Since A — XI : 
Dom(A) — > H is one-one and onto, A is an eigenvalue of A + 7CZ) if and only if 
— 1/7 is an eigenvalue of CD(A-XI)^ 1 . Both implications in the first sentence now 
depend on the elementary fact that if U, V are vector spaces over C, X : U — > V, 
Y : V — >■ U are linear operators and a G C is non-zero, then a is an eigenvalue of 
XY if and only if it is an eigenvalue of YX. □ 



In spite of the second statement in Lemma 20, the £(/C)-valued function m is 
easier to analyze than the scalar function p. One says that the analytic function 
m : C + —> C{K) is an operator-valued Herglotz function if (m(X)f,f) G C + for 
every / G /C\{0} and A G C+; see [3]. 
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Lemma 21 Suppose that B = B* > 0, C = D = B 1 / 2 , K. is the closure of the 
range of B and * is the operation of truncation to /C. Then 

m (A) = (B 1/2 (A - Xiy 1 B 1/2 ) ii (10) 
is a £(lC)-valued Herglotz function. 

Proof The assumptions imply that g = B 1 ! 2 f and that 

(m(X)fJ) = ((A-XI)- 1 g,g}, 

which lies in C + by the spectral theorem. □ 

Theorem 22 If B = B* > has finite rank N and A G C + then there are at least 
1 and at most N values of ^ such that X is an eigenvalue of A^; all such 7 lie in 

c+. 

Proof If / e Dom(A), / ^ and A y f = Xf then 

(AfJ) +1 (BfJ}=X(fJ}. 

This implies that 

Im (7) {Bf, f) = Im(A) (/, /) . 
Since the right hand side is positive, we deduce that (Bf, f) > and Im(7) > 0. 



Lemma 20 states that A is an eigenvalue of A^ if and only if — 1/7 is an eigenvalue 



of the N x N matrix m(A). The final statement of Lemma 21 implies that every 
eigenvalue of m(A) lies in C+. This proves that there are at least 1 and at most iV 
distinct values of 7, each of which lies in C+. □ 



The equation (11) below is a special case of the Nevanlinna-Riesz-Herglotz repre- 



sentation of operator- valued Herglotz functions; see [3j. 



Lemma 23 Under the assumptions of Lemma 21, let P(E) denote the spectral 



projection of A associated with any Borel subset £CR, If 

Q(E) = (B 1 / 2 P{E)B l ' 2 f 

then Q is a finite, non-negative, countably additive, £(JC)-valued measure on R 
and 

m ( X )= f ^tW s ) (u) 



R 3 ~ X 



for all X G C + . 
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Proof The formula (11) follows directly from 



JR s - A 

which is proved using the spectral theorem. □ 
The following lemma is well-known, but we include a proof for completeness. 

Lemma 24 If f : [a, b] — » C is bounded and measurable then 

r-b 



f(s)Q(ds)\\<\\f\\ 00 \\Q([a,b})\\. 

Proof If </>, if) e /C, we have 

|( I* f(s)Q(ds)<t>M = |( f b f(s)P(ds)(B 1 / 2 <l>),(B 1 ^))\ 



= \(f(A)P([a,b))(B^<f>),(B^)}\ 

= |(/(A)(P([a,6])i? 1 /^ );( p (K6])i? i/^ )) | 

< ||/|| 00 ||P(K6]) J B 1 /^|||| P([a)6])jB i^|| 

= ||/|| oo (g([a ) 6])0 ) 0) 1/2 (g([a,6])^^) 1/2 

< ll/||oo||g(k6])||||0||||V||. 

The lemma follows. □ 



Our next lemma defines two operators that will be used in Theorem 26 



Lemma 25 Let [a, b] C R. Then there exist bounded, self-adjoint operators X, Y 
K, — > K such that 



X = / Q(ds), (12) 

J a 

X 1/2 YX 1/2 = / sQ(ds). (13) 



a 



Moreover < X < and al <Y < bl . 

Proof Since X = Q([a, b}), the inequalities for X and the boundedness of X follow 
from 

< <Q([o, &])/,/) < (Q(R)fJ) = (B*f,f) < \\B*\\ " "' 2 
valid for all / G JC. If one defines 

Z= I sQ(ds) 

J a 
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then aX < Z < bX. If X is invertible this immediately implies that 

al < Y = X- 1/2 ZX- 1/2 < bl. 

The general case follows as in Lemma [2] □ 

Determining the eigenvalues of m(A) for a given range of values of A can sometimes 
be aided by writing 

m(A) = rrii(X) + 7712(A) 

where mi may be computed more readily than m and /77.2(A) can be neglected or 
replaced by an appropriate approximation for the selected range of values of A. 
Theorem 26 enables one to replace the contribution of an interval [a, b] to 777(A) by 
a single operator provided A is far enough away from [a, b] . 

Theorem 26 Let 

m ( X )= I -^Q(ds), 



where X G C + and Q is a finite, non-negative, countably additive, C(K.)-valued 
measure on R. Given [a, b] C R and L > 0, let 

772(A) = / ^g(d s )+X 1 / 2 (F-A/)- 1 X 1 / 2 , (14) 

Js<£{a,b] S — X 

where X and Y are as defined in Lemma 25_, Then 

|m(A)-m(A)|<^^||Q([a,fc])|| 
for all X G C + such that dist(A, [a, b}) > L. 
Proof It suffices to prove that 

II [ j^~ x Q(ds) - X^(Y - XI)-W\\ < 2 -^j^\\Q([a,b])\\ (15) 

for all A satisfying the stated conditions. 
We first observe that 

1 1 s-b (s-bf 



s-X X-b (X-b) 2 (s-X)(X-b) 2 ' 
Integrating both sides with respect to Q over [a, b] yields 

" b 1 X X^YX^-bX 



(16) 



r 1 x 



(A - b) 2 

6 (s - b) 2 

Q(ds), (17) 



a (s-X)(X-b) 2 
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and then 



s-X 



< 



X-b 
(b-a) 2 



(a - by 



L 3 



IIQ(MDII 



by Lemma 24 



We next use the formula 
I 



/ Y-bl 

+ 1 7 + 



(y - biy 



to obtain 



Y — XI X-b (X-b) 2 (Y-XI)(X-b) 2 



\XV 2 (Y - XI)- 1 X 1 ' 2 + + X 1 ' 2 * b L x^\ 

X — o 



(19) 



< 



(x-by 

ii X i/2 {y-bi? 1/2|l 

11 (F-A/)(A-6) 2 11 
(6 - a) 2 



L 3 



-||Q(k 



The proof of (15) is completed by combining (18) and (20). 



(20) 
□ 



Remark 27 One can obtain a better approximation than that in Theorem 26 
if [a, b] is divided into several subintervals, each of which is used to produce an 



extra term in the formula (14). The new rh is, of course, more cumbersome to use 
numerically. □ 

From this point onwards we assume that A is a possibly unbounded self-adjoint 
operator and that v4 7 = A + where Bf = (f, e)e and e G H is a vector with 



norm 1. The Herglotz function (10) is then scalar- valued and given by the formula 
m(A) = ((A - XI)~ l e } e). 

If one approximates m uniformly in a chosen region by another analytic function 
whose zeros are more easily computed, then one can apply Rouche's theorem to 
approximate the zeros of m and hence the spectrum of A 1 . Lemma 28 is directly 
applicable to the polynomial p(j, X) = det(v4 7 — XI), where 7 is fixed. The connec- 



tion between this and the Herglotz function m is explained in Lemma 15 and (24) 



The proof of Lemma 28 can be adapted to cases in which p is not a polynomial; one 
needs an upper bound on the orders of its zeros, a lower bound on the distances 
between zeros and a lower bound on p for points that are not close to a zero. 

Lemma 28 Let < e < 1/2, let U be a bounded open set in C and let U £ = {z G 
C : dist(z,U) < 2e}. Let p be a monic polynomial such that \p(z)\ > e for all 
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z G U £ \U . Suppose that every root of p is simple and that \s — 1\ > 2 for any two 
distinct roots of p. Let q be an analytic function on U £ such that \p(z) — q(z)\ < e 
for all z G U £ . If A G U £ and p(\) = then A G U and there exists exactly one zero 
of p and one zero of q inside the circle C\ j£ = {w G U £ : \w — A| = s}. If X G U £ 
and q(\) = then A G U and there exists exactly one zero of p and one zero of q 
inside C\ >£ . 

Proof The assumptions of the lemma imply immediately that neither p nor q can 
vanish in U £ \U. Suppose that A G U and p(A) = 0. Then \s — A| > 2 for all s in 
the finite set S of roots of p that are not equal to A. If w G C\ :2e then \w — s\ > 1 
for all s G S. Therefore 



\p{w)\ = 



.w 



s&S 



W 



> \w - A| = 2s. 



Since C\^ £ and its interior are contained in U £ , we may apply Rouche's theorem 
to p and q on and inside Ca,2 £ , and deduce that q has exactly one zero inside C\^ £ - 
The same argument evidently applies to C\ t£ . 

On the other hand if A G U and q{\) = then |p(A)| < s. If T is the set of all 
roots of p then 



£>|p(A)| 



teT 



so I A — 1\ < 1 for at least one root of p; from this point onwards we use the symbol 
t to refer to one such root. If S — T\{t} then |A — s\ > 1 for all s G S, so t is 
unique. Therefore 



e>|p(A)| = 



(A-t)H(A- S ) 



ses 



>\X-t\ 



Repeating the first paragraph of the proof with A replaced by t, both p and q have 
exactly one root inside the circle C t) 2e- This contains the region inside C\ t£ , so both 
p and q have at most one root inside C\ t£ . The proof is concluded by noting that 
we have already observed that they have at least one root inside C\ >£ . □ 

Example 29 Given positive integers M 1; M 2 and L, let iV = Mi + M 2 and let A 
be the diagonal N x N matrix with entries 



l Mi+s,Mi+s 



= (r-l)/(M!-l), 
= L + 1 + (s - 1)/(M 2 



1), 



for 1 < r < Mi and 1 < s < M 2 , so that Spec(A) C [0, 1] U [L + 1, L + 2]. Also 
let Ay = A + 7_B where B is the rank one operator associated with the unit vector 
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e E C N denned by e r = (2M 1 )~ 1 / 2 for 1 < r < M 1 and e Ml+s = (^M 2 )- X l 2 for 
1 < s < M 2 . One sees immediately that ||e|| = 1 in C^. 

The continuous curves in Figure [2] show parts of six of the spectral curves of 
for 7 = te ie when M l = 5, M 2 = 25, L = 4, < t < 20 and = 89°. Most of the 
curves stay within a small distance of their starting point as t increases. The curve 
starting at the eigenvalue 0.5 of the 30 x 30 matrix A moves rapidly away from the 
real axis as t increases but eventually converges to 3. There is only one curve that 
diverges to oo as t — > oo, and a part of this appears in the top right-hand part of 
the figure. 

The dashed curves in Figure [2] are produced in a similar manner but with Mi = 5 



and M 2 = 1, so that A is a 6 x 6 matrix. Following the prescription of Theorem 26 
we define A rr as above for 1 < r < 5, but put A 66 = 5.5, so that the function 

(10) 



m of Theorem 
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is the Herglotz function for the pair A, e, where e r = (lO) -1 / 2 
for 1 < r < 5 and e§ = 2 -1 / 2 . In spite of the substantial reduction in the size of 
the matrix, the part of the spectrum in {A : Re(A) < 2} is almost unchanged, as 
predicted by Theorem 26 □ 




Figure 2: Spectral curves in Example 29 
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6 Rank one perturbations 



In this section we obtain more detailed results of the type already considered under 
the assumptions that A is a self-adjoint operator acting on the Hilbert space H 
and that Bf = (f, e)e for all / G H, where e is a unit vector in %. 

We define A 7 on Ti by 

A,f = Af + 1 (f,e)e (21) 

where 7 G C. We summarize a few of the many results known in the case 7 G R 
and then consider non-real 7, for which new issues arise. We also assume that e is 
a cyclic vector for A in the sense that ||e|| = 1 and C = lm{A n e : n = 0, 1, . . .} is 
dense in "H. This is equivalent to B being cyclic for A by Corollary [7j 

The four propositions below provide the general context within which our more 
detailed results are proved. The first is classical and may be found in [8]. 

Proposition 30 Let e be a cyclic vector for the bounded self-adjoint operator A 



and let v4 7 be defined by (21). J/ 7 G R then every eigenvalue of has multiplicity 
one. If a G R and X a is an isolated eigenvalue of A a , then X a can be analytically 
continued to all real 7 that are close enough to a and A 7 > for all such 7. 

Proposition 31 If A G C\Spec(A) and 7 G C then A is an eigenvalue of A^ if 
and only if 

l+7((A-A/)- 1 e,e> = 0. (22) 

This formula defines 7 as an analytic function of A G C\Spec(A); one has 7(A) = 
— l/m(A), where 

m(X) — {(A — A/) _1 e, e). (23) 



This is a special case of Lemma 20 In finite dimensions one may alternatively use 
the formula 

= l + 7 ((A-A/)- 1 e,e). (24) 

See Lemma dU 



Proposition 32 T/ie function m(A) defined for all A G C\Spec(A) fry (S3) is 
a Herglotz function in the sense that m(A) G C± /or a// A G C±. Moreover 
\m(x + iy)\ < 1/M /or a// x G R and ?/ ^ 0. If A is bounded then m(A) 7^ for 
all A G C sfic/i t/ia£ |A| > ||A||. It follows that 

7(A) = A+(Ae,e)+0(|A|- 1 ) 

as |A| —7- 00. 
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Proposition 33 Suppose that 70, Ao G C + satisfy (22) and that Xq has algebraic 
multiplicity 1 as an eigenvalue of A 10 . Then there exists an analytic function X 
0/7 defined for all 7 in some neighbourhood 0/70 such that 7 and A satisfy (22). 
Moreover A' (7) ^ in this neighbourhood. 

Proof The first statement of the proposition is a standard fact from perturbation 
theory for the eigenvalues of operators that depend analytically on a parameter. 
Given this, we differentiate (|22j) with respect to 7 to obtain 

((A - XI)-' e, e) - 7 A'( 7 ) ((A - XI)~ 2 e, e) = 0. 

Assuming that the neighbourhood is small enough, 7 ^ R and ((A— XI)~ 1 e, e) 7^ 
by Proposition 32 This implies that A' (7) 7^ 0. 

□ 

In the rest of this section we assume that iV = dim('H) < 00, that e G H has norm 
one and is a cyclic vector for A, and that Im(7) > 0. Our goal is to understand 
how the eigenvalues of A^ depend on 7 for very small and very large 7, and the 
mapping properties from the one asymptotic regime to the other. We start with 
the case in which 7 is real and positive. 

Lemma 34 Under the assumptions of the last paragraph, let Ai,...,Ajv be the 
eigenvalues of A written in increasing order and let Sx, ■ ■ ■ , 5v-i be the eigenvalues 
of the truncation A^ of A to K, = {/ : (/, e) = 0}. Then 



Ai < 5i < X 2 < ■ ■ ■ < 5n-i < A 



A 1 - 



If one assumes that 7 > then the eigenvalues X rn of A^ are all strictly increasing 
analytic functions 0/7 satisfying A rj0 = A r . Moreover lim^+oo A n7 = 5 r for 1 < 
r < N — 1 and lim 7 _;. +00 Ajv, 7 = +00. 



Proof This uses Proposition [30] and the variational formula for the eigenvalues of 

Ay. □ 

We now turn to the study of the case 7 G C+. 

Theorem 35 Given 9 G (0, it), define 

S e = |JSpec(A tel9 ). (25) 
t>o 

Then 

S e n S$ = if 6 £ <f) (26) 

and 

U S e = C+. (27) 

0e(O,7r) 

Moreover the limit set of each Se in C + U {00} is Spec(A) U Spec^) U {00}. 
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Proof If A G C + then (22) determines 7 = te 10 uniquely. This fact implies (26) and 
(27). The limit set of Sq is the union of lim^o Spec(Aj e ie) and limt_>. +o0 Spec(A te ie) , 
both of which were determined in Theorem [TH □ 



We now turn to the structure of the individual sets Sg. Let Si, ... , <5at_i denote the 
eigenvalues of A\ written in increasing order and let 6~n = 00. As before we say 
that a curve a : (0, 00) — > C + is simple and analytic if it is a one-one, real analytic 
mapping and a'{t) is non-zero for all t G (0, 00). 

Theorem 36 There exists a finite increasing set T C (0, it) such that if £ 
(0,7r)\T then Sg is the union of N disjoint simple analytic curves. Each curve 
starts at some A r G Spec (A) and ends at some 8 t m where r is a permutation of 
{1, 2, . . . , iV}. This permutation is constant in each subinterval J of (0, ir)\T, but 
it may change from one interval to another. 



Proof This a corollary of Theorems [19] and [35j but some of the calculations are 
simpler because the polynomial p(7, A) defined in ^ has the following explicit 
form. By expanding the determinant using an orthonormal basis whose first term 
is e one obtains 



p(7,A) = det(A- A/)+7det(^ 
= Po(A)+tpi(A), 



•A/*) 



(28) 
(29) 



where ^ denotes the truncation to K = {(f) : (0, e) 
degree N and p\ is a polynomial with degree iV - 



0}, po is a polynomial with 



1. The formula (29) can also 



be derived from (22). Let F be the finite exceptional set defined just before The- 
orem 19 It follows from (26) and (27) that there is a finite set T C (0,7r) such 
that 9 G T if and only if te*^G F for some t > 0. If # G (0, 7r)\T then the curve 
g(t) = te %e lies in Qq, as defined just before Theorem 19, which yields most of the 



statements of this theorem. 6 ^ T implies that the curves are simple and non- 
intersecting because every A G C + is associated with only one 7 = te td and hence 
with only one value of t G (0, 00) by (22). The constancy of the permutation on 



each subinterval J follow from the continuous dependence of the curves in Sg on 9. 
To prove the last statement, it is sufficient to consider the following example. □ 



Example 37 Let H = C 2 and let 



A, 



1 +7« 2 
70J/3 



7«/3 
-1 + 7/3 2 



where 7GC, a>0, (3 > and a 2 + (3 2 = 1, so that e = 
is a cyclic vector for A = Aq. The eigenvalues of Ay are 



(a, P) has norm one and 



A 



7 



±,7 



± 



T 



+ l- 7 (/3 2 -a 2 ). 
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Figure [3] plots these eigenvalues for a = a/3/2, /3 = 1/2, 7 = te %e and < t < 4. 
The two dashed curves correspond to the choice 9 = 119°; the two intersecting 
continuous curves correspond to the choice 9 = 120°; the two dotted curves corre- 
spond to the choice 9 = 121°. Note that the only critical point is 7 C 
so T = {120°}. The corresponding eigenvalue of A lc is A c = { — 1/2 + za/3/2}, 
which has algebraic multiplicity 2 but geometric multiplicity 1, by a direct compu- 



tation or Theorem 10 It is clear that the permutation r of the set {1, 2} defined in 
Theorem [36] is different for 9 < 120° and for 9 > 120° and that there is no natural 
way of defining such a permutation for 9 = 120°. □ 




Figure 3: Spectral curves described in Example 37 



Example 38 Let A be the 5x5 diagonal matrix with eigenvalues X r = r, r = 
1, 2, 3, 4, 5, and let e = (1, 1, 1, 1, l)/\/5- Then the limits of the eigenvalues of A 7 
as I7I 00 are given numerically by fii = 1.35556, /12 = 2.45608, /i3 = 3.54390, 
//4 = 4.64442 and /is = 00. For each 9 exactly one of the five eigenvalue curves 
diverges to 00. The table below lists the permutations r associated with each angle 
9 E (0°, 180°) that is a multiple of 10°. 
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e 


r(l) 


r(2) 


r(3) 


r(4) 


r(5) 


0° 


1 


2 


3 


4 


5 


10° 


1 


2 


3 


4 


5 


20° 


1 


2 


3 


4 


5 


30° 


1 


2 


3 


4 


5 


40° 


1 


2 


3 


4 


5 


50° 


1 


2 


3 


4 


5 


60° 


1 


2 


3 


4 


5 


70° 


1 


2 


3 


5 


4 


80° 


1 


2 


3 


5 


4 


90° 


1 


2 


5 


3 


4 



The first change in r occurs for Q\ G (61°, 62°), while the second occurs for 6 2 G 
(81°, 82°). □ 



7 The limit N -> oo 

The previous analysis clarifies to some extent how the spectra of a family of N x N 
matrices A 1 = A+^B depend on 7 for very small and very large 7. However, it does 
not capture the full range of phenomena that can occur for 7 of intermediate sizes. 
Even if one is interested in a particular fairly large value of N, one often obtains 
further insights by considering a family of N x N matrices = An + 7-B/v- 
From this point of view the case N = 00 is regarded as an idealization that may 



be simpler to analyze than the original problem. Results such as Proposition 40 
may then be used to estimate the difference between the two cases. 

We assume throughout that B N f = (/, ejv)ejv for all / G where G is a 
unit vector. The set of eigenvalues of A Na is obtained by solving 1 + •ym N (X) = 0, 
or equivalently 7 = — l/mjv(A), where m N are the Herglotz functions 

m N (\) = ((A N - AIjv) _1 ejv,ejv). 



See Propositions 31 and 32 



Theorem 39 Let 

Hn{i) = max {Im(A r) Ar >7 ) : 1 < r < N} 

where 7 G C + and {A rj Ar i7 }^ 1 are the eigenvalues of A Nn repeated according to 
their algebraic multiplicities. Then 

^(7) > Im( 7 )/iV (30) 
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for all N > 1 and 7 G C+. Suppose that \\An\\ < c /or a// iV and i/iai converge 
to moo locally uniformly on C + as iV — » 00. IfjE C + and Ra.n(m 00 ) i/ien 

lim ^(7) = 0. 

N—toc 



Proof We note that is a Herglotz function, unless it is a real constant. Each 
function : C + — > C + is surjective because mjv(A) = — 1/7 has N solutions 
A G C + counting multiplicities, namely the eigenvalues of A^ a . We show in 



Example 41 that need not be surjective. The lower bound (30) follows from 



AW7) > Im X r,N,-y J = Im (tr(Ajv + 7#a0) = Im( 7 ). 

If 7 G C + then every eigenvalue A n jv i7 of v4tv, 7 satisfies 

|A r ,iv )7 | < \\A N \\ + |t| ll-Bivll < c + |7|. (31) 

Therefore < Hn(i) < c + |t|- If A* at (7) does not converge to as iV — > 00 then 
there exists a subsequence iV(s) and a constant C2 > such that ^n{s)(j) > C2 for all 
s; and then a subsequence r(s) such that Im(A r ( s ) j jv( s ), 7 ) > C2 for all s. Combining 
this with (31), there exist subsubsequences, which we again parametrize using s, 
and A G C + such that lim^^ X r ( s ),N(a),'y — X where Im(A) > c 2 and |A| < c + |7|. 
Since m N ^(X r ^ N ^^) = — 1/7 for all s, the local uniform convergence of m N to 
moo implies that moo (A) = —1/7. □ 

Theorem [39] depends on the assumption that m n converges to N 00. 

The following proposition allows one to estimate the difference between mjv(A) and 
m 00 (A) for problems of the above type by putting 

|2 



k(s) 



\9W 



A 



where the choice of g depends on the problem. It may be seen that the bound on 
the difference is 0(Im(A)~ 2 ) as Im(A) — > 0. 

Proposition 40 Let k be a continuous function on [a, b] with bounded first deriva- 
tive and let N be a positive integer. Then 



u N 
— a 



N 

where 



J2k(a + r(b-a)/N) = ! k(s) ds + 

T = l •' Cl 



rem (32) 



rem < \\k 

1 2N 
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Proof The left hand side of (32) is the sum of N terms of the form 

•cr d 



5k(c r 



c r —5 



ds 



[(s — c r + 6)k(s)] ds 



k(s) ds + C r 



c r —5 



where 5 = {b — a)/N, c r = a + r5 and 



lev I 



(s — c r + 5)k'(s)) ds 



< \\k'\ 



S\ds 



\k'\ 



Summing over r one obtains 



rem < \\k 



{b-af 
2N 



\k'\ 



□ 



Example 41 Let An be the N x N diagonal matrix with entries -Ajv.n.n = n/N for 
all n and let the rank one matrix associated with the unit vector = N~ l l 2 
for all n. Then Spec^^r) C [0, 1] and 



m N 



^ N^n/N-X' 

r=l ' 



It may be seen that converges locally uniformly to 



moo(A) 



1 ds 



log 



o s — X \ X 



A - 1 



(33) 



as iV — > oo. It follows that the range of moo is {z G C : < Im(z) < ir} and the 
set of fi G C + such that moo (A) = — I// 2 has no solution is the closed disc 

D = { 7 :| T -i/(27r)|<l/(27r)}. 



Figure |4] provides a contour plot /^jv(t) for N = 100, the contours corresponding 
to the values 0.01, 0.02, 0.03, 0.04, 0.05 of ^(7). The circle dD is included for 
comparison. □ 
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Figure 4: Contour plot of ^n(i) for N = 100. 



Example 42 Let be the N x N diagonal matrix with entries A^^ n ^ n = n/N 
for all n and let the rank one matrix associated with the unit vector ejv,n = 
N~ l l 2 g{n/N) where g G L 2 (0, 1) and g is sufficiently regular. Then Spec(AAr) C 
[0, 1] and rriN converges locally uniformly to 

™-(A) = /' (34) 

7n s - A 



as N — > oo, where f(s) = \g(s)\ . The integral in (34) is well-defined for every 
A G C + because / G -^ 1 (0, 1), but the form of the range of depends on whether 
either or both of the integrals 

1 — S 



is finite. In Example 41 both integrals are infinite. More generally the first integral 
diverges if and only if A + 'jB has a negative eigenvalue for all real negative 7, 
while the second integral diverges if and only if A + ^B has a positive eigenvalue 
for all real positive 7 

The range of moo is the union of the sets m 00 (S r ejT .), where 

S etr = {A : Im(A) > e and |A| < r} . 



30 



These increase monotonically as s > decreases to and as r increases to oo. The 
boundary dS £>r may be parametrized as a simple closed curve 7 £jr . A standard 
theorem in complex analysis states that each set m 00 (S' e . r ) is the union of the 
range of the closed curve a £ ^ r = o 7 e>r and the set of all z not in this range 
whose winding number with respect to cr £jr is non-zero. 

The observations above allow one to compute the range of moo approximately by 
taking e > small enough and r large enough. This is particularly easy if one can 
write moo in closed form. If g(s) = s 1 ^ 2 for all s G [0, 1] then 

™oo(A) = / S . ds = 1 + A log 
Jo s - x 

Figure [Hjwas obtained by putting e = 1CT 8 . The set of 7 for which — 1/7 = moo(A) 
is not soluble is the part of C+ that is inside the closed curve plotted. This curve 
starts at —1 and ends at 0. The gap observed near is a numerical artifact that 
arises because the convergence to is logarithmic. □ 



0.75 - 



0.5 - 




Figure 5: Boundary curve in Example [42] 
Theorem 43 Suppose that 

m (A) = [ -f^-ds (35) 
Jk s - A 
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for all A G C +; where f(s) > for all seR, 



, ds < oo 

and < H/Hoo < oo. Tnen m(C + ) is contained in 

{z:0<lm(z) <7r||/|| 0O }. 
Hence —1/7 = m(A) is not soluble for any 7 in £/ie closed disc 

1 



7 : 



7 



2^ 



< 



2tt 



(36) 



(37) 



(3? 



Proof The condition (36) ensures that the integral (35) defining m(A) converges 



for all A G C + and defines an analytic function of A. The condition f(s) > and 
< H/Hoo ensures that the range of m is contained in C+. 

We next observe that if A = u + iv where u G R and v > then 



Im(m(A)) 



R 



(n — s) 2 + 



:/(s)ds. 



A direct estimate of this integral yields (|37|) , and (38) follows. 



□ 
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